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Abstract 



In previous papers we showed that the data of the DAMA/Nal experi- 
ment for WIMP direct detection, which indicate a possible annual modulation 
effect, are widely compatible with an interpretation in terms of a relic neu- 
tralino as the major component of dark matter in the Universe. In the present 
note we discuss the detectability of the relevant supersymmetric neutralino 
configurations by two of the most promising methods of indirect search for 
relic particles: measurement of cosmic-ray antiprotons and measurement of 
neutrino fluxes from Earth and Sun. 
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I. INTRODUCTION 



In Refs. |T]-§| we showed that the indication of a possible annual modulation effect in 
WIMP direct search [J|,|5] are interpretable in terms of a relic neutralino which may make 
up the major part of dark matter in the Universe. 

We recall that the DAMA/Nal data H single out a very delimited 2-<r C.L. region in 
the plane £^caiar 0n ' ) ~ m x' where m x is the WIMP mass, CTg™^ 011 ^ is the WIMP-nucleon 
scalar elastic cross section and £ = p x /pi is the fractional amount of local WIMP density p x 
with respect to the total local dark matter density p\. In Fig. 1 we display the region R, 
which is obtained from the original 2-cr C.L. region of Ref. || (where p\ is normalized to the 
value pi = 0.3 GeV cm -3 ), by accounting for the uncertainty in the value of pi. In obtaining 
region R we have conservatively considered the range: 0.1 GeV cm -3 < pi < 0.7 GeV cm -3 , 
because of a possible flattening of the dark matter halo and a possibly sizeable baryonic 
contribution to the galactic dark matter ||. In Fig. 1 is also displayed the scatter plot of 
the supersymmetric configurations which, in our scanning of the susy parameter space, turn 
out to be contained in region R (the set comprised of these configurations is called set S). 
The theoretical framework adopted here is the Minimal Supersymmetric extension of the 
Standard Model (MSSM), as illustrated in Refs. to which we refer for details. 

In |]-f3| we have proved that a significant fraction of configurations of set 5* is explorable 
at accelerators in the near future. In the present paper, we discuss the detectability of 
these configurations by different means, i.e. by two of the most promising ways of indirect 
search for relic particles: measurement of cosmic-ray antiprotons [0HTT| and measurement 



of neutrino fluxes from Earth and Sun 12-14 



Preliminary discussions on the links between the possible annual modulation effect and 
the results of indirect searches have been presented in Refs. [0,0], in connection with the 
former DAMA/Nal data (data collected in the running period # 1) @- Here we extend our 
investigation to the new results of the DAMA/Nal experiment (running period #2) |5|], which 
establish a further indication for the annual modulation effect and, combined with the for- 
mer ones, single out an effect at 99.6% C.L. In the present paper we also discuss neutralino 
cosmological properties in terms of its relic abundance Q x h 2 . This quantity is evaluated 
here as illustrated in [ [L5[| , and for the parameter £ we take £ = min(f2 x /i 2 /(f2/i 2 ) min , 1), 
with (f2/i 2 ) min = 0.01. In exploring our supersymmetric parameter space we apply an up- 
per bound, conservatively set at the value £l x h 2 < 0.7, and we consider as cosmologically 
interesting the range 0.01 < Q x h 2 < 0.7. However, we stress that, according to the most 
recent data and analyses fl6| , the most appealing interval for the neutralino relic abundance 
is 0.02 < tt Y h 2 < 0.2. 

Most of our results for the various quantities considered in this paper will be provided in 
the form of scatter plots, obtained by varying the representative points of the susy parameter 
space over the configurations of set 5*. When in a specific figure (or in a panel of it) the value 
of p\ is fixed, it is meant that the corresponding scatter plot concerns only the subset of S 
comprised of those configurations which belong to that particular value of pi. For instance, 
for pi = 0.3 GeV cm -3 , only the points inside the dashed curve of Fig. 1 are considered. 
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II. COSMIC RAY ANTIPROTONS 



Relic neutralinos, if present in our halo, would produce antiprotons by pair annihilation 
{LT| . To discriminate this potential source of primary cosmic-ray p's from the secondary 
antiprotons, i.e. p's created by interactions of primary cosmic-ray nuclei with the interstellar 
medium, one can use the different features of their low-energy spectra (T p < 1 GeV, T p being 
the antiproton kinetic-energy). In fact, in this energy regime the secondary p spectrum drops 
off very markedly because of kinematical reasons, while antiprotons due to neutralinos would 
show a milder fall off. 



In Ref. Ill] the low-energy cosmic-ray p spectrum has been calculated, in terms of a 
possible contribution of cosmic-ray antiprotons due to neutralino annihilation in the galactic 
halo and of a newly calculated flux for secondary antiprotons. The total spectrum has 
been compared with the data collected by the balloon-borne BESS spectrometer during 
its flight in 1995 (hereafter referred to as BESS95 data) Hl7j . The BESS95 data consist 



of 43 antiprotons grouped in 5 narrow energy windows over the total kinetic-energy range 
0.180 GeV < T p < 1.4 GeV. 

In the present paper we evaluate the cosmic-ray antiproton flux as in Ref. []rTJ , restricting 
the supersymmetric configurations to those of set 5* and we compare our predictions with 
the BESS95 data. We refer to [TT|] for all the details pertaining the evaluation of the 
production of primary and secondary antiprotons as well as for the properties related to 
their propagation in the halo and in the heliosphere. Here we only recall the features of the 
neutralino mass distribution function adopted in JIT] as well as here. This mass distribution 



function is taken spheroidal and parameterized as a function p x (r, z) of the radial distance r 
from the galactic center in the galactic plane and of the vertical distance z from the galactic 
plane in the form 

Px<r,z)= Px Q2+ % + + r ° 2//2 , (1) 

where a is the core radius of the halo, r Q is the distance of the Sun from the galactic center 
and / is a parameter which describes the flattening of the halo. Here we take the values: 
a = 3.5 kpc, r Q = 8 kpc. The quantity p x denotes the local value of the neutralino matter 
density, which is factorized as p x = £pi, in terms of the total local dark matter density p\. 
The density p\ is calculated by taking into account the contribution given by the matter 
density of Eq.(0) to the local rotational velocity. We use here the value pi = 0.3 GeV cm -3 
in the case of a spherical halo (/ = 1). When / < 1 (oblate spheroidal distribution), pi is 
given by ||0 

p,(f) = p l (f=l) JJ ^= Ji . (2, 

Our results will be shown for the following values of the local density: pz/(GeVcm -3 ) = 
0.1 (/ = 1), 0.3 (/ = 1), 0.5 (/ = 0.50), 0.7 (/ = 0.33). For each value of p h we will give the 
top-of-atmosphere (TOA) antiproton fluxes, as the sum of the secondary flux and of the 
primary flux due to neutralino annihilation for the various supersymmetric configurations 
of set S, pertaining to that specific value of pi. 
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Fig. 2 displays the scatter plots for the top-of-atmosphere (TOA) antiproton fluxes 
calculated at Tp = 0.24 GeV, to conform to the energy range of the first bin of the BESS95 
data (0.175 GeV < Tp < 0.3 GeV). This figure also shows the band of the experimental 
results in this bin. We find that, while most of the susy configurations of the appropriate 
subset of S stay inside the experimental band for pi = 0.1, 0.3 GeV cm -3 , at higher values 
of pi a large number of configurations provide p fluxes in excess of the experimental results. 
This occurrence is easily understood on the basis of the different dependence on p\ of the 
direct detection rate and of the p flux, linear in the first case and quadratic in the second 
one. 

To present the comparison of our theoretical fluxes with the experimental ones in a more 
comprehensive way, we give in Fig. 3 the scatter plots of the reduced chi-square Xred values, 
obtained for the configurations of each of the four subsets of S, in a fit to the BESS data 
over the whole energy range (0.175 GeV < Tp < 1.4 GeV). In these scatter plots we have 
included only configurations which are at least at the level of the experimental value (within 
1-er) in the first energy bin. We remark that, apart from values of p\ close to the largest 
value of its physical range, a large sample of the susy configurations, relevant for annual 
modulation data in direct search of relic neutralinos, provides a good fit to the current data 
on cosmic antiprotons. 

By way of example, we show in Fig. 4 a sample of fluxes obtained with a few susy 
configurations of set S which give a good fit to the BESS95 data. The parameter P, which 
denotes the various supersymmetric configurations, is the gaugino fractional weight (for 
instance, a neutralino which is a pure gaugino has P = 1). 

From Fig. 3 it also turns out that, especially at large pi, many configurations of set S have 
very large values of x^ ed . Configurations which have xLi > 4 may be considered as strongly 
disfavoured on the basis of the present experimental data on cosmic antiprotons (xLi = 4 
corresponds to a 99.9% C.L. in our case of 5 degrees of freedom). However, because of the 
many uncertainties involved in the current theoretical evaluations (for instance, effects due 
to solar modulation of the original interstellar fluxes) and in the experimental measurements 
of the antiproton spectrum, a strict exclusion criterion may be applied only to configurations 
which are in large excess of x^d = 4- By the same arguments, from our results we derive that 
high values of pi, i.e. pi ~ 0.7 GeV cm -3 , appear to be disfavoured. Hereafter, configurations 
with x;?cd > 4 will be denoted by special symbols. 

In Fig. 5 we give x? c d versus the neutralino relic abundance. On the basis of these 
results, it turns out that the central values of the range for the local density, 0.3 GeV cm -3 
< pi < 0.5 GeV cm -3 , is the best suited for neutralinos of cosmological interest. The break at 
Q x h 2 = 0.01 in the shape of the scatter plots in this figure is due to the rescaling in the local 
density introduced through the parameter £ for Q x h 2 < (Qh 2 ) m i n = 0.01. Because of this, 
one roughly has: $p oc (o" ann ) _1 oc Q x h 2 for VL x h 2 < (Qh 2 ) m i n , and $p oc <r ann oc (f2 x /i 2 ) -1 
otherwise, where a ann is the neutralino-pair annihilation cross section. 

The various results discussed in this section show the remarkable property that a number 
of the supersymmetric configurations singled out by the annual modulation data may indeed 
produce measurable effects in the low-energy part of the p spectrum. Thus the joint use 
of the annual modulation data in direct detection and of the measurements of cosmic-ray 
antiprotons is extremely useful in pinning down a number of important properties of relic 
neutralinos and show the character of complementarity of these two classes of experimental 
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searches for particle dark matter. This stresses the great interest for the analyses now under 
way of new antiproton data, those collected by a recent balloon flight carried out by the 
BESS Collaboration and those measured by the AMS experiment [^5j during the June 
1998 Shuttle flight. 



III. NEUTRINOS FROM EARTH AND SUN 

The signals we are going to discuss now consist of the fluxes of up-going muons in a 
neutrino telescope, generated by neutrinos which are produced by pair annihilations of neu- 
tralinos captured and accumulated inside the Earth and the Sun [12|-[14||. The process goes 



through the following steps: a) capture by the celestial body of the relic neutralinos through 
a slow-down process due essentially to neutralino elastic scattering off the nuclei of the 
macroscopic body, b) accumulation of the captured neutralinos in the central part of the 
celestial body, c) neutralino-neutralino annihilation with emission of neutrinos, d) propaga- 
tion of neutrinos and conversion of their component into muons in the rock surrounding 
the detector (or, much less efficiently, inside the detector), and finally e) propagation and 
detection of the ensuing up-going muons in the detector. 

The various quantities relevant for the previous steps are calculated here according to 
the method described in Ref. |T3| , [14|j , to which we refer for further details. Here we only wish 



to recall that the neutrino fluxes are proportional to the annihilation rate of the neutralinos 
inside the macroscopic body, Ta, which is given by [|T^] 

c . 2 rt 



-„ 2 tanh^-j, (3) 

where t is the age of the macroscopic body (t ~ 4.5 Gyr for Sun and Earth) and ta = 
(CCa)^ 1 ^ 2 ', Ca is the annihilation rate proportional to the neutralino-neutralino annihilation 
cross-section and C is the capture rate, proportional to the neutralino-nucleus cross-section 
and to the neutralino local density. From Eq. (||) it follows that in a given macroscopic body 
the equilibrium between capture and annihilation (i.e. Ta ~ C/2 ) is established only when 

t > T A - 

From the evaluation of the annihilation rate for neutralinos inside the Earth and the 
Sun it turns out that, for the Earth, the equilibrium condition depends critically on the 
values of the model parameters, and is not realized in wide regions of the parameter space. 
Consequently, for these regions the signal due to neutralino annihilation may be significantly 
attenuated. Capture by the Earth is very efficient for m x < 70 GeV because of mass- 
matching condition between the neutralino mass and the nuclear mass of some important 
chemical constituents of the Earth (O, Si, Mg, Fe). On the contrary, in the case of the Sun, 
the capture-annihilation equilibrium is reached for the whole range of m x , due to the much 
more efficient capture rate due to the stronger gravitational field [pOfl . 

Let us now report our results. In Fig. 6 we display the scatter plots for the flux of the 
up-going muons from the center of the Earth, for various values of the local dark matter 
density pi. In this figure is also reported the current 90% C.L. experimental upper bound on 
$Ei arth , OD t a i nec i by MACRO [^1|] (a similar experimental upper bound is given in Ref. |[22|| ). 
We note the particular enhancement at m x ~ (50 — 60) GeV, due to the mass-matching 
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effect. It is also noticeable an effect of suppression and spreading of the fluxes for m x > 70 
GeV at pi > 0.3 GeV cm -3 . This is due to the fact that the configurations with these 
large values of p\ may imply, because of the annual modulation data, a neutralino-nucleus 
cross-section which is too small to establish an efficient capture rate, necessary for the 
capture-annihilation equilibrium in Earth. This is shown in Fig. 7, whose scatter plots show 
that indeed in the case of large values of p\ many configurations induce sizeable deviations 
from the value tanh 2 (t/t^) = 1, which guarantees equilibrium in the macroscopic body. 

In Fig. 6 the configurations which would be excluded on the basis of the antiproton data, 
as previously discussed in Sect. II, are denoted differently from those which would survive 
this criterion. 

By comparing our scatter plots with the experimental MACRO upper limit, one notices 
that a number of supersymmetric configurations provide a flux in excess of this experimen- 
tal bound and might then be considered as excluded. However, it has to be recalled that 
a possible neutrino oscillation effect []23| may be operative here and affect the indirect neu- 
tralino signal as well as the background consisting in atmospheric neutrinos. The influence 
of a neutrino oscillation on the problem under discussion would depend very critically on 
the nature of the oscillation itself (v^ converting either into a v T or into a sterile neutrino). 
Therefore a strict enforcement of the current upper bound on $^ arth should be applied with 
caution as long as the neutrino oscillation properties are not fully understood. However, 
it is rewarding that the set 5* of supersymmetric configurations is quite accessible to relic 
neutralino indirect search by measurements of up-going fluxes, as displayed in our results 
of Fig. 6. 

Fig. 8 displays the up-going muon fluxes from the Sun. We notice that a number of 
configurations of higher masses are easily accessible to the measurements of these fluxes. 

IV. CONCLUSIONS 

We have proved that a sizeable fraction of the supersymmetric neutralino configurations 
singled out by the DAMA/Nal data on a possible annual modulation effect in WIMP direct 
search, may provide signals detectable by two of the most interesting methods of indirect 
search for relic particles: measurement of cosmic-ray antiprotons and measurement of neu- 
trino fluxes from Earth and Sun. 

For the case of cosmic antiprotons, it has been shown that present data are well fitted 
by total spectra which include a p contribution from neutralino-pair annihilation, with 
neutralino configurations which are relevant for annual modulation in direct detection. These 
data can also be used to reduce the total sample of the supersymmetric configurations 
under study, and to narrow the range of the local density, by disfavouring its largest values. 
Investigation by measurements of cosmic p's looks very promising in view of the collections 
and analyses of more statistically significant sets of data in the low-energy regime which are 
currently under way |24|,|25| and which may soon provide further relevant information. 

Measurements of neutrino fluxes from Earth and Sun, due to capture and annihilation 
of neutralinos inside these celestial bodies, have been proved to be sensitive to neutralino 
configurations singled out by the annual modulation data. An appropriate interpretation 
of these measurements preliminarily requires some clarification of the oscillation neutrino 



6 



properties. 

In conclusion, we have shown how the measurements of cosmic-ray antiprotons and 
of neutrino fluxes from Earth and Sun are complementary to those of direct detection in 
investigating astrophysical and cosmological properties of the relic neutralinos which may 
be at the origin of the possible modulation effect in direct search. 
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FIGURE CAPTIONS 



Figure 1 - Scatter plot of set S in the plane ^ x _ £ a scaiar • The dashed contour line 
delimits the 2-er C.L. region, obtained by the DAMA/Nal Collaboration, by combining 
together the data of the two running periods of the annual modulation experiment ||. The 
solid contour line is obtained from the dashed line, which refers to the value pi = 0.3 GeV 
cm -3 , by accounting for the uncertainty range of pi, as explained in Sect. I (the region 
delimited by the solid line is denoted as region R in the text). Displayed in this figure are 
only the representative points of the susy parameter space, which fall inside the region R (this 
set of configurations is defined as set S). Dots and circles denote neutralino configurations 
depending on their values for the neutralino relic abundance. 

Figure 2 - Scatter plots for the TOA antiproton fluxes calculated at Tp = 0.24 GeV 
versus the neutralino mass. The configurations of set S are subdivided into the 4 panels, 
depending on the corresponding value of the local density: p//(GeVcm~ 3 ) = 0.1 (/ = 1), 
0.3 (/ = 1), 0.5 (/ = 0.50), 0.7 (/ = 0.33). The dashed lines denote the central value and 
the 1-cr band of the BESS95 data in the first energy bin: 0.175 GeV < T p < 0.3 GeV. 

Figure 3 - Scatter plots of the Xred values, obtained for the configurations of each of 
the four subsets of S, in a fit to the BESS95 data over the whole energy range (0.175 GeV 

< Tp < 1.4 GeV) versus the neutralino mass. Only configurations which are at least at 
the level of the experimental value (within 1-cr) in the first energy bin are included. The 
three horizontal lines represent: the Xred i n a ^ of the data using the secondary flux only 
(dot-dashed line), and the Xred i n a fit using the total (primary plus secondary) flux at 95% 
C.L. (short-dashed line) and at 99.9% C.L. (long-dashed line). The pi values are in units 
of GeV cm~ 3 . 

Figure 4 - Some examples of TOA fluxes which provide a good fit to the BESS95 data, 
for pi = 0.3 GeV cm -3 . The dotted line denotes the secondary antiproton flux only. The 
solid, dashed and dot-dashed lines denote the total antiproton fluxes, obtained as a sum 
of the secondary spectrum and of the primary spectra due, respectively, to the following 
neutralino configurations: m x = 62 GeV, Q x h 2 = 0.03 and P = 0.98 (solid line); m x = 81 
GeV, Vl x h 2 = 0.02 and P = 0.99 (dashed line); m x = 95 GeV, Vt x h 2 = 0.03 and P = 0.99 
(dot-dashed line). 

Figure 5 - Scatter plots of the Xred values, obtained for the configurations of each of 
the four subsets of S, in a fit to the BESS data over the whole energy range (0.175 GeV 

< Tp < 1.4 GeV) versus the neutralino relic abundance. Notations as in Fig. 3. 

Figure 6 Scatter plots for the up-going muon fluxes from the center of the Earth versus 
the neutralino mass. The configurations of set S are subdivided into the 4 panels, depending 
on the corresponding value the local density: 

P;/(GeVcm" 3 ) = 0.1 (/ = 1), 0.3 (/ = 1), 0.5 (/ = 0.50), 0.7 (/ = 0.33). Dots denote 
configurations which could be excluded on the basis of the BESS95 antiproton data, circles 
denote configurations which survive this exclusion criterion. The dashed line denotes the 
MACRO upper bound. 

Figure 7 - Scatter plots of the quantity tanh 2 (t / ta) versus the neutralino mass. Other 
specifications are as in Fig. 6. 

Figure 8 - Scatter plots for the up-going muon fluxes from the Sun versus the neutralino 
mass. Other specifications are as in Fig. 6. 
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